function plot_gould_flux5m(t1,t2)
% plot_gould_flux5m(t1,t2) plots the 5-min
% wind stress and heat flux data stored in
% gould_met5m for the period [t1:t2]. 
% t1,t2 are optional input. default is 
% t1=min(t), t2=max(t).
% 4/11/01

% get data in t1-t2
load gould_met5m
if nargin==0
   t1=min(t);t2=max(t);
end
if nargin==1
   t2=t1+1;
end
m=find(t>=t1 & t<=t2);
t=t(m);stress=stress(m);qnet=qnet(m);qsw=qsw(m);
qlw=qlw(m);qsen=qsen(m);qlat=qlat(m);wd=wd(m);
d2r=pi./180;
tx=stress.*cos(d2r*wd);
ty=stress.*sin(d2r*wd);

% setup figure
clf
orient tall
t0=floor(min(t));t1=ceil(max(t));

subplot(6,1,1)
plvec(t,tx+i*ty)
ylabel('Tau')
axis([t0 t1 -1.5 1.5]);
title('R/V Gould 01-03 5-min Surface Fluxes')

subplot(6,1,2)
plot(t,stress,'b')
ylabel('Tau')
grid
axis([t0 t1 0 1.5]);

subplot(6,1,3)
plot(t,wd,'b')
ylabel('Dir')
grid
axis([t0 t1 -200 200]);

subplot(6,1,4)
plot(t,qnet)
ylabel('Qnet')
grid
axis([t0 t1 -500 500]);

subplot(6,1,5)
plot(t,qsw,'b',t,qlw,'r')
ylabel('Qsw, Qlw(red)')
grid
axis([t0 t1 -200 700])

subplot(6,1,6)
plot(t,qsen,'b',t,qlat,'r')
ylabel('Qsen, Qlat(red)')
grid
axis([t0 t1 -200 50]);

xlabel('2001 Yearday')
hold off